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ABSTRACT 

We  present  a  modified  beat-frequency  modulated  accretion  model  appropriate 
to  a  description  of  the  horizontal  branch  oscillations  of  Z-type  low-mass  x-ray 
binaries.  This  model  incorporates  several  new  physical  elements:  (i)  magnetic 
screening  by  the  accretion  disk  and  neutron  star,  (ii)  radiation  pressure 
dominance  in  the  inner  disk,  and  (iii)  stellar  rotation  effects.  Fitted  to 
observational  data,  the  model  constrains  the  surface  magnetic  field  in  a  typical 
LMXB  Z  source  to  be  less  than  10®  G  and  the  stellar  spin  frequency  to  be 
greater  than  400  Hz.  The  predicted  high  stellar  spin  frequency  falls  in  a  range 
where  observational  upper  limits  on  pulsation  are  not  a  constraint.  The  derived 
magnetic  field  and  spin  frequency  range  are  compatible  with  the  fastest  radio 
pulsars. 
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1.  Introduction 

Quasiperiodic  oscillations  (QPO)  with  intensity  dependence  (van  der  Klis  et  al.  1985) 
are  a  common  feature  in  Z-sources.  These  sources  are  low  mass  x-ray  binaries  (LMXBs) 
whose  x-ray  emission  evolves  through  a  characteristic  sequence  of  spectral  ratios,  as 
observed  in  a  color-color  diagram.  The  three  spectral  modes  comprising  the  “Z”  are  known 
as  the  horizontal,  normal,  and  flaring  branches.  In  these  sources  the  correlation  of  QPO 
frequency  with  source  brightness  manifests  itself  when  the  source  is  on  the  horizontal 
branch;  for  this  reason  the  oscillations  are  also  called  “horizontal  branch  oscillations” 
(HBOs).  This  effect  distinguishes  them  from  other  QPO  seen  when  the  source  is  in  one  of 
the  other  spectral  states. 

A  model  for  HBO  has  been  proposed  (Alpar  and  Shaham  1985)  in  which  the  QPO 
frequency,  i/q,  is  the  beat  frequency  of  the  Kepler  orbit  frequency,  z/^-,  at  the  magnetosphere 
radius  and  the  neutron  star  spin  frequency,  i/^,  i.e., 

Uq  =  VK{Rm)  -  I's  (1) 

A  more  detailed  model  by  Lamb  et  al.  (1985)  specifies  how  luminosity  oscillations  at  vq  are 
produced  through  modulation  of  the  mass  transfer  rate  to  the  neutron  star  in  a  transition 
zone  at  the  inner  disk  boundary.  The  papers  cited  above  lay  out  the  essential  features  of 
the  beat  frequency  modulated  accretion  (BFMA)  model. 

A  fit  of  parameters  from  standard  magnetosphere  models  gives  a  magnetic  field 
~  6  X  10®  G  (or  magnetic  moment  /z  =  6  x  lO^^G  cm^)  and  a  spin  frequency  z/,  ~  110 
Hz  in  the  case  of  GX  5-1.  Similar  values  of  B  and  can  be  inferred  for  other  sources 
from  QPO  data  obtained  while  the  source  is  on  the  horizontal  branch.  The  magnetospheric 
modeling  concepts  used  in  the  BFMA  model  derive  from  the  model  of  Ghosh  and  Lamb 
(1979;  hereafter  GL)  for  accreting  binary  pulsars.  However,  the  GL  magnetospheric  model 
originally  pertained  to  high  mass  x-ray  binaries  (HMXB)  with  large  magnetospheres.  For 
LMXBs,  it  needs  modification  for  effects  arising  when  the  magnetosphere  radius  is  near 
the  neutron  star  surface  and  radiation  pressure  modifies  the  disk  structure.  In  addition, 
relativistic  effects  in  rapidly-spinning  neutron  stars  that  can  affect  angular  momentum 
transfer  in  the  system  need  to  be  considered.  There  are  also  technical  fine  points.  For 
example,  in  the  GL  torque  theory  the  fastness  parameter,  =  i^sli^K{Bm),  plays  a  central 
role.  For  a  certain  critical  value  of  this  parameter,  a;,c,  the  accretion  torque  vanishes.  The 
numerical  value  of  Wjc  used  in  binary  pulsar  models  is  not  necessarily  valid  in  LMXBs  and 
has  even  been  questioned  in  high-mass  binary  pulsars  (Wang  1987,  Shibazaki  1989).  We 
present  here  an  alternative  BFMA  model  which  addresses  these  problems.  Our  results 
suggest  that  the  spin  periods  in  LMXB  Z  sources  are  much  shorter  than  previously  thought, 
and  the  magnetic  moments  are  significantly  smaller. 
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Summaries  of  the  successes  of  the  BFMA  model  and  also  of  some  issues  concerning  it 
may  be  found  in  Lamb  (1989)  and  van  der  Klis  (1989).  The  former  reference  shows  how  the 
BFMA  models  for  HBO  can  be  unified  with  a  broader  treatment  of  all  behavior  of  Z-type 
LMXB  sources.  The  latter  reference  is  a  more  comprehensive  review  of  QPO  observations. 
The  BFMA  successfully  fits  the  intensity  dependence  of  the  HBO.  The  transition  zone 
modulation  picture  accounts  for  the  correlation  of  the  QPO  intensity  with  low-frequency 
noise  (LFN),  although  there  remains  an  unresolved  issue  concerning  how  this  correlation 
extends  to  the  shortest  timescales  (Norris  et  al.  1990,  Mitsuda  et  al.  1991).  The  reader  is 
also  referred  to  Lewin,  Van  Paradijs  k  Van  der  Klis  (1988)  for  a  review  of  QPO  observations 
in  a  number  of  sources  and  a  description  of  other  models  of  the  QPO  phenomena. 

A  major  success  for  the  beat  frequency  model  is  an  observation  of  a  QPO  in  the 
transient  source  EXO  2030+375  (Angelini,  Stella  and  Parmar  1989)  which  conforms 
phenomenologically  to  the  beat  frequency  model  with  parameters  re-scaled  into  a  different 
range  of  z/,  and  B  appropriate  to  accreting  binary  pulsars.  Because  this  source  is  detected 
as  a  pulsar,  the  spin  period  is  a  known  quantity  in  this  case  and  not  "a  free  parameter  to  be 
fitted  to  observations  such  as  variation  of  QPO  frequency  with  intensity.  This  case  argues 
persuasively  that  intensity  modulation  at  the  beat  frequency  between  VK{Rm)  and  i/,  is 
a  reality  in  some  neutron  star  binary  systems.  Validation  of  the  standard  BFMA  model 
is  less  complete  in  the  LMXB  systems  for  which  it  was  first  proposed.  Indeed  there  are 
significant  problems  with  it  in  that  context.  Most  conspicuously,  careful  searches  have  not 
detected  coherent  x-ray  pulsations  in  those  systems  (Wood  et  al.  1991,  Hertz  et  al.  1990, 
Vaughan  et  al.  1994).  Upper  limits  on  coherent  modulation  near  or  below  1  percent  have 
been  established  for  spin  frequencies  up  to  500  Hz  in  all  6  of  the  Z-type  sources  (see  Table 
1).  Our  goal  is  to  achieve  a  satisfactory  match  of  the  BFMA  picture  to  observations  of 
HBO  in  LMXBs. 

An  additional  motivation  to  re-examine  the  estimated  spin  frequencies  of  these  systems 
is  the  possibility  of  an  evolutionary  link  from  low-mass  x-ray  binary  systems  to  radio 
millisecond  pulsairs  (Helfand,  Ruderman  k  Shaham  1983,  Joss  k  Rappaport  1983,  Paczynski 
1983,  Savonije  1983).  The  fastest  of  the  millisecond  pulsars  have  spin  periods  near  1.6  ms 
and  magnetic  fields  as  low  as  10*  G  (Camilo,  Thorsett  k  Kulkarni  1994),  values  significantly 
lower  than  those  obtained  by  conventional  BFMA  model  fits  to  HBO.  We  will  explore 
whether  such  radio  pulsars  are  compatible  with  our  modified  BFMA  picture. 

The  spin  period  of  LMXB  neutron  stars  is  a  central  issue.  Therefore,  we  start  with  a 
more  complete  discussion  of  the  interpretation  of  upper  limits  obtained  from  period  searches. 
Subsequent  sections  present  the  modified  BFMA  model,  apply  it  to  HBO  observations, 
discuss  alternative  ways  to  treat  details,  and  finally  return  to  further  observational  issues 
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and  tests. 


2.  Interpretation  of  Upper  Limits  on  Pulsation 

Wood  et  al.  (1991)  summarize  why  LMXB  spin  periods  should  lie  in  the  millisecond 
range.  The  case  for  this  conclusion  was  persuasive  even  before  QPO  were  discovered.  The 
picture  of  disk  accretion  with  at  most  a  weak  magnetic  field  requires  eventual  spin-up  to 
high  frequency.  Accreting  material  continually  adds  angular  momentum  to  the  star,  except 
perhaps  near  breakup  (Popham  &  Narayan  1991).  Detection  of  HBO  in  a  source  makes 
the  case  even  stronger  if  the  BFMA  model  is  accepted.  However,  sensitive  searches  for  spin 
periods  in  several  LMXB  QPO  sources  have  not  detected  coherent  modulation  down  to  the 
1%  level  in  pulsed  amplitude  (Vaughan  et  al.  1994). 

The  interpretation  of  limits  on  pulsations  in  Z  sources  has  been  presented  in  several 
earlier  papers,  primarily  in  the  context  of  determining  whether  limits  for  spin  frequencies 
near  100  Hz  were  consistent  with  the  standard  BFMA  model  or  possibly  constituted  a 
basis  for  rejecting  that  model.  We  now  re-examine  the  issue  from  a  different  standpoint, 
taking  the  BFMA  picture  as  provisionally  correct  and  asking  instead  what  is  the  most 
reasonable  value  for  the  spin  frequency.  We  describe  several  new  points  not  covered  in  the 
older  discussions. 

Pulsation  upper  limits  have  been  growing  more  stringent  with  refinements  in  search 
technique  and  improvements  in  data.  When  the  BFMA  was  first  proposed  the  upper  limits 
for  some  Z  sources  were  near  10%.  More  recent  upper  limits  reach  levels  of  a  few  times 
10~^  out  to  frequencies  approaching  500  Hz  (see  Table  1),  which  is  the  Nyquist  frequency 
for  observations  conducted  with  the  Ginga  satellite  (Wood  et  al.  1991,Vaughan  et  al.  1994). 
At  higher  spin  frequencies  it  also  becomes  progressively  more  expensive  computationally  to 
overcome  smearing  of  the  Fourier  power  spectrum  due  to  the  unknown  orbital  acceleration, 
which  makes  searching  trial  grids  of  acceleration  parameters  a  necessity. 

Several  mechanisms  have  been  proposed  to  account  for  the  non-detections  of  spin 
periods.  Since  a  magnetic  field  and  channeling  of  flow  near  the  star  are  central  to  the 
BFMA  picture  it  is  difficult  to  avoid  uneven  heating  of  the  surface  that  results  in  production 
of  emission  modulated  at  the  spin  frequency.  The  mechanisms  proposed  to  reduce  this 
pulsation  below  observable  levels  include  (i)  gravitational  lensing  and,  (ii)  scattering  in  a 
hot  cloud  or  corona  surrounding  the  star.  We  review  each  of  these  in  turn. 

Gravitational  lensing.  Wood  et  al.  (1988)  calculated  the  effects  of  lensing,  exploring 
variation  of  polar  cap  sizes,  surface  gravity,  and  geometry  of  the  observer  and  field  axis. 
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Two  antipodal  caps  were  used,  a  choice  consistent  with  a  dipole  field  and  one  that  minimizes 
resulting  pulsation.  A  neutron  star  with  a  radius  of  2  Rs  (8.3  km  for  IAMq)  can  reduce 
pulsation  to  about  7  x  10“^  relative  amplitude  modulation  for  a  45°  misalignment  between 
the  dipole  field  axis  and  the  observer’s  line-of-sight.  Such  a  small  radius  is  consistent  with 
only  the  softest  of  the  proposed  equations  of  state  for  neutron  star  matter  (Friedman, 
Ipser  &  Parker  1986).  If  the  stellar  radius  is  increased  to  2.5/?*  (10.4  km)  then  the  pulse 
amplitude  is  reduced  only  to  about  1.4%. 

These  numbers  are  consistent  with  the  limits  available  in  1987  but  more  recent  limits 
reach  lower  values,  hence  the  lensing  explanation  is  becoming  marginal  even  for  slowly 
rotating  stars.  A  further  problem  has  been  raised  by  Chen  and  Shaham  (1989),  who  have 
studied  effects  of  extremely  high  surface  velocity  (rapid  rotation).  Most  of  their  analysis 
pertains  to  very  high  velocities,  such  as  would  be  found  near  spin  periods  of  1  millisecond. 
Special  relativistic  beaming  effects  and  redshifting  from  the  emission  frame  to  the  observer’s 
frame  acting  on  realistic  emission  spectra  can  overcome  the  gravitational  lensing  effects  and 
enhance  pulsation.  The  difficulty  in  applying  these  results  to  neutron  stars  with  parameters 
implied  by  the  standard  BFMA  model  is  that  the  surface  velocities  are  very  small  compared 
with  those  used  in  the  calculations,  only  0.02  c  at  the  surface  of  a  star  with  radius  10  km 
and  a  spin  period  of  10  ms.  The  effects  increeise  roughly  eis  the  surface  velocity. 

The  primary  uncertainty  in  evaluating  the  effects  of  gravitational  lensing  on  pulse 
supression  comes  from  our  ignorance  of  the  quantity  M/R  for  neutron  stars.  Available 
theoretical  models  of  neutron  stars  are  not  very  constraining  and  observational  evidence  is 
also  inconclusive.  If  future  observational  results  should  settle  this  point  in  favor  of  stellar 
radii  greater  than  about  9  km  (for  1.4M®)  then  gravitational  lensing  as  an  explanation  for 
the  lack  of  detectable  coherent  pulsations  will  be  excluded,  leaving  only  the  explanations 
that  follow. 

Coronal  scattering.  Bussard  et  al.  (1988)  investigated  the  effects  of  a  hot,  spherical 
scattering  cloud  on  QPO  behavior.  They  found  that  power  spectra  intrinsic  to  the  source 
can  be  severely  attenuated  at  high  frequencies  by  scattering  in  a  hot  coronal  cloud  with 
the  amount  of  attenuation  depending  on  optical  depth  and  cloud  radius.  The  intrinsic 
power  spectra  can  be  related  to  the  observed  spectra  by  a  simple  frequency- dependent 
transmission  function  T{v)  =  Fobs  I  Pinitiai  of  the  form 

TM  =  e-^'(l  +  +  4xV=)]  (2) 

At  high  frequencies,  i/  >  (7c,/27r)e’',  this  function  asymptotically  approaches  where 
r  is  the  Thompson  optical  depth.  At  lower  frequencies  the  transmission  function  has  a 
roughly  Lorentzian  shape,  falling  off  as  above  a  frequency  corresponding  to  'yal2T, 
where  7^  is  approximately  given  by  the  inverse  of  the  product  of  the  optical  depth  of  the 
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cloud  and  its  light-crossing  time.  Early  spectral  fittings  (White  &  Holt  1982,  Shibazaki  & 
Mitsuda  1984,  ,  Hirano  et  al.  1984,  )  suggested  a  hot  cloud  with  optical  depth  ranging  from 
less  than  1  to  a  few.  More  recent  fittings  of  Exosat  data  for  4  Z-sources  estimate  optical 
depths  in  the  range  4-7  (Schulz  &  Wijers  1993).  Observed  temporal  lags  in  the  QPO  source 
Cygnus  X-2  (Hasinger  1986,  Mitsuda  &  Dotani  1989),  if  due  primarily  to  scattering,  suggest 
a  cloud  radius  greater  than  1000  km  (Bussard  et  al.  1988). 

As  pointed  out  by  Bussard  et  al.,  their  results  strictly  apply  only  to  an  isotropically 
emitting  source  which  pulses  in  time.  Coherent  pulsations  from  a  QPO  source  would  likely 
arise  due  to  beaming  from  the  surface  of  a  rapidly-spinning  neutron  star  in  which  case 
both  isotropization  of  the  beam  due  to  angular  diffusion  and  smearing  due  to  variations  in 
photon  residence  time  in  the  cloud  must  be  considered.  Brainerd  and  Lamb  (1987)  treated 
the  former  effect.  However  if  the  stellar  rotation  period  is  sufficiently  short  compared  to  the 
diffusion  time,  tRc/c  ~  7“^  then  radial  diffusion  will  be  the  most  important  effect  and  the 
Bussard  et  al.  results  can  be  used  to  estimate  the  attenuation  of  coherent  pulsations  (as 
well  as  QPO  oscillations)  by  the  cloud. 

In  the  standard  BFMA  model  (Alpar  and  Shaham  1985,  Lamb  et  al.  1985),  the  neutron 
star  rotation  frequencies  are  estimated  to  be  just  above  100  Hz  and  the  QPO  frequency 
ranges  from  5  to  as  high  as  50  Hz.  As  shown  by  Bussard  et  al.  (1988),  if  we  require  that 
the  power  at  the  QPO  frequency  be  relatively  undiminished  while  power  above  100  Hz  be 
strongly  attenuated,  the  radius  of  the  scattering  cloud  is  constrained.  That  is  uq  <  (7a/27r). 
For  r  =  5  this  implies  Rc  <  550  km.  At  frequencies  above  (7a/2x)  the  power  attenuation 
is  approximately  (7a/27r)^/i/^  until  it  reaches  a  limiting  value  of  at  high  frequencies. 
Thus  if  UQmax  ~  50  Hz,  then  the  coherent  pulse  amplitude  is  reduced  to  0.50  its  original 
value  for  i/g  =  100  Hz  and  0.08  for  1/3  =  640  Hz. 

In  Figure  1  we  compare  the  pulsational  upper  limits  for  LMXB  sources  with  the 
attenuation  expected  from  coronal  scattering  and  from  gravitational  lensing.  The  viability 
of  either  of  these  modulation  supression  mechanisms  is  beginning  to  be  severely  strained  by 
the  observed  upper  limits  for  frequencies  below  approximately  500  Hz.  It  may  be  that  the 
simplest  explanation  for  the  lack  of  detection  of  spin  periods  in  LMXB  sources  is  that  the 
spin  periods  are  less  than  2  ms.  In  the  remainder  of  this  paper  we  explore  the  viability  of 
this  hypothesis.  Ultimately,  however,  the  issue  will  be  resolved  observationally. 
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3.  Magnetosphere  Model 

The  relationship  giving  the  magnetosphere  radius  in  terms  of  the  mass  accretion  rate 
M  and  the  stellar  magnetic  moment  fj,  in  the  conventional  BFMA  model  is 

Rm  =  0.52Ra  =  34km  (3) 

where  /X27  and  Mn  are  the  stellar  magnetic  moment  and  mass  accretion  rate,  respectively, 
in  units  of  lO^^G  cm^  and  lO^^gm  s~^.  m  is  the  neutron  star’s  mass  in  units  of  lAf0,  and 
Ra  is  the  spherical  Alfven  radius. 

Although  there  is  minor  disagreement  as  to  the  value  of  the  numerical  prefactor,  most 
authors  assume  this  dependence  (Alpar  and  Shaham  1985,  van  den  Heuvel  1977,  Lamb 
et  al.  1985).  The  relation  was  applied  by  GL  to  magnetically-coupled  accretion  torques 
acting  on  x-ray  pulsars  in  binary  systems.  Magnetic  moments  in  those  cases  are  typically 
much  larger  than  lO^^G  cm^.  Spin  periods  are  on  the  order  of  the  Kepler  period  at  R^-, 
seconds  to  hundreds  of  seconds.  In  LMXB  systems  magnetic  moments  are  weaker  by 
193-4  'j'jjjg  causes  the  magnetosphere  to  be  far  smaller,  hence  new  physical  conditions  and 
considerations  arise. 

The  resemblance  of  HBOs  in  LMXBs  to  the  one  seen  in  the  HMXB  binary  pulsar 
EX02030-I-375  shows  it  is  reasonable  to  attempt  to  extend  these  ideas  to  the  smaller 
LMXB  magnetospheres.  That  is  the  framework  for  our  model,  i.e.,  we  assume  that  a  disk, 
a  magnetosphere,  and  a  relatively  narrow  disk-magnetosphere  transition  zone  continue  to 
exist  in  the  LMXB  case.  However,  at  least  three  factors  cause  assumptions  carried  forward 
from  the  HMXB  disk-magnetosphere  picture  to  become  questionable  in  this  limit,  (i)  Spin 
frequencies  become  very  high  (i.e.  Vj  approaches  UK{Rm))  and  can  even  approach  the 
breakup  or  mass-shedding  limit,  (ii)  Because  the  disk  likely  comes  within  a  few  stellar 
radii  of  the  neutron  star’s  surface,  magnetic  screening  can  distort  the  field  configuration 
substantially  from  a  dipole,  (iii)  When  the  accretion  luminosity  Lx  approaches  the 
Eddington  limit  the  accretion  disk  just  outside  the  magnetosphere  is  likely  dominated  by 
radiation  pressure,  ultimately  becoming  unstable  when  the  luminosity  approaches  very  close 
to  or  exceeds  Lsdd-  However,  our  working  hypothesis  is  that  for  LMXB  Z  sources  on  the 
HB  the  luminosity  is  low  enough  that  the  required  revision  is  an  incorporation  of  radiation 
pressure  into  the  disk-magnetosphere  model  and  not  some  more  radical  recourse. 

Some  of  these  needs  for  revision  have  been  noted  in  earlier  work.  For  example. 
White  and  Stella  (1988)  have  offered  prescriptions  for  including  radiation  pressure  in  the 
determination  of  the  magnetosphere  radius  within  the  context  of  the  model  of  GL.  They 
consider  a  number  of  different  models  and  argue  that  none  can  be  excluded  as  models  for 
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at  least  some  of  the  observed  QPOs.  However,  not  all  these  models  are  workable  for  HBO 
because  not  all  give  intensity  dependence,  to  say  nothing  of  the  other  phenomenology. 

Aside  from  the  revisions  to  the  conventional  picture  of  the  disk-magnetosphere 
interaction  that  are  necessary  for  bright  LMXB  sources,  there  are  other,  more  generic, 
deficiencies  in  the  conventional  magnetosphere  models.  For  example,  Wang  (1987)  hcis 
pointed  out  an  inconsistency  in  the  treatment  of  the  magnetic  pitch  angle  in  the  disk  in 
the  original  Ghosh-Lamb  model.  Wang’s  revision  of  the  GL  model  leads  to  a  substantially 
larger  value  of  the  critical  fastness  parameter  u>sc,  primarily  because  Wang  predicts  a  smaller 
value  for  the  pitch  angle  in  the  outer  disk  and  therefore  a  smaller  spin-down  contribution  to 
the  total  accretion  torque  on  the  star.  If  the  QPOs  observed  in  HMXB  sources  are  due  to 
the  beat  frequency  mechanism,  then  the  fastness  parameters  derived  from  the  observations 
are  in  reasonable  agreement  with  Wang’s  model  applied  to  gas  pressure  dominated  disk 
accretion  (see  discussion  in  Appendix).  However,  Wang’s  model  does  not  have  a  narrow 
transition  zone  between  Keplerian  flow  in  the  disk  and  the  approximately  co-rotating 
magnetosphere.  This  zone  is  a  feature  of  the  Ghosh-Lamb  model  and  appears  to  be  a 
necessary  feature  if  the  beat  frequency  model  of  QPO  is  correct. 

There  are  a  number  of  other  problems  with  any  of  the  conventional  disk-magnetosphere 
models.  Given  the  theoretical  difficulties  of  the  problem  and  the  present  state  of 
observations,  a  compelling  theory  that  is  applicable  to  LMXBs  is  not  likely  to  be  available 
soon.  Nonetheless,  it  is'  possible  to  make  progress  on  the  spin  period-QPO  dilemma  in 
LMXBs  by  the  following  approach:  (i)  we  introduce  a  set  of  modifications  to  the  standard 
treatment  of  disk-magnetosphere  interactions  to  deal  with  known  deficiencies;  (ii)  we 
parameterize  these  modifications  and  combine  them  into  a  “modified  BFMA  model”,  and 
(iii)  we  explore  this  parametrized  model  over  a  range  of  the  parameters  that  we  believe 
is  reasonable,  to  see  how  various  modifications  constrain  the  spin  frequency  and  magnetic 
moment  of  neutron  stars  in  Z  sources.  This  is  done  by  requiring  agreement  of  eq.(l)  with 
the  intensity-dependent  QPO  frequency  observed  on  the  horizontal  branch.  Anticipating 
the  conclusion,  we  find  that  over  a  substantial  part  of  the  plausible  range  of  the  parameters 
of  the  model,  the  spin  frequencies  inferred  are  much  larger  than  those  implied  by  the 
conventional  model  and  that  the  neutron  star  magnetic  moments  required  are  about  the 
same  as  those  of  the  fastest  millisecond  radio  pulsars.  Additionally,  the  new  model  offers 
some  insights  into  issues  regarding  the  QPO  shot  rate  and  the  correlation  of  QPO  power  and 
low  frequency  noise  (LFN)  power  raised  by  Ginga  observations  (Norris  et  al.  1990,  Mitsuda 
et  al.  1991).  While  we  believe  these  modifications  are  physically  well-motivated,  it  should 
be  noted  that  at  this  time  it  is  virtually  impossible  to  directly  establish  the  validity  of  this 
modified  model. 
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The  details  of  the  modifications  to  the  disk-magnetosphere  interaction  picture  are 
presented  in  the  remaining  subsections.  Here  we  only  summarize  these  modifications  and 
the  associated  parameters.  The  reader  not  interested  in  the  details  can  skip  to  section  4. 
The  areas  of  modification  and  the  associated  parameters  are: 

(i)  Angular  Momentum  Transfer:  Throughout  the  remainder  of  this  section  the  disk  is 
modeled  with  radiation  pressure  as  the  dominant  pressure  in  the  inner  regions,  according 
to  the  prescription  summarized  by  Treves  et  al.  (1988).  Torques  at  the  magnetosphere 
boundary  are  modeled  keeping  terms  of  order  S/Rm,  where  6  is  the  width  of  the  transition 
zone.  Vertical  mass  loss  from  the  transition  zone,  associated  with  material  transfered  onto 
closed  field  lines,  is  explicitly  taken  into  account  with  a  parameter  /dm,  which  is  unity  when 
there  is  no  vertical  loss. 

(ii)  the  magnetic  pitch:  We  consider  the  original  model  of  GL  in  which  the  magnetic 

pitch  angle  70  =  5^/  Bz  in  the  narrow  transition  zone  is  of  order  unity.  We  also  consider 
the  model  of  Wang  (1987)  in  which  the  pitch  angle  depends  on  physical  parameters  of  the 
disk  and  the  star.  We  parameterize  7^  in  our  model  by  70  =  (7,^)^^,  where  0  <  <  1  and 

7u,  is  Wang’s  expression  for  the  pitch  angle  given  by  eq.  (8)  below. 

(iii)  width  of  transition  layer:  We  essentially  follow  the  Ghosh-Lamb  treatment.  An 

important  parameter  in  determining  the  width  of  the  transition  zone  is  the  average  radial 
velocity  Vr  through  the  transition  zone.  GL  take  Vr  to  be  the  free-fall  velocity  u/y,  given  by 
eq.  (14)  below.  However,  when  the  fastness  parameter  approaches  unity,  Vr  may  be  better 
approximated  by  the  radial  drift  velocity  in  the  disk  Vd  just  outside  the  transition  zone.  Vd 
depends  on  the  viscosity  of  the  disk.  In  order  to  explore  a  range  of  models  between  these 
limits  we  parametrize  the  radial  velocity  as  Vr  =  (u//)^~^'’(u<i)^‘',  where  0  <  <  1. 

(iv)  effective  magnetic  moment:  We  model  the  bare  stellar  magnetic  field  as  a  dipole 
with  moment  p.  As  originally  noted  by  GL,  screening  currents  in  the  disk  and  the 
diamagnetic  behavior  of  the  neutron  star  will  lead  to  an  enhancement  of  the  field  in  the 
magnetosphere.  While  this  is  probably  negligible  for  x-ray  pulsars,  in  the  case  of  LMXBs, 
where  p  is  likely  to  be  several  orders  of  magnitude  smaller  and  the  disk  extends  to  within  a 
few  stellar  radii  of  the  surface,  the  enhancement  effect  is  more  important.  We  model  this 
enhancement  in  terms  of  an  effective  dipole  moment  pgff{r).  A  parameter  of  importance 
that  appears  later  is  A^  =  —{Rm/peff)dpeff/dRm- 
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3.1.  Angular  Momentum  Transfer 


Following  the  approach  of  GL,  we  derive  an  expression  for  the  magnetosphere  radius 
that  explicitly  includes  the  effects  of  stellar  rotation.  We  assume  that  the  Keplerian  flow 
of  the  disk  ends  at  the  radius  where  the  magnetic  coupling  to  the  disk  is  large  enough  to 
remove  sufficient  angular  momentum  from  the  flow  to  enforce  co-rotation.  The  relevant 
equation  of  angular  momentum  conservation  is 


d{Mr^Q,{r)) 


=  r  U— - h  -|-  viscous  stresses. 

dv 


Next,  we  assume  that  the  transition  from  nearly  Keplerian  flow  in  the  disk  to  co-rotation 
occurs  in  a  narrow  transition  boundary  layer  of  width  S  «  r.  Ignoring  viscous  stresses 
(since  magnetic  stress  >>  viscous  stress  in  the  transition  layer),  we  have,  after  integrating 
and  keeping  terms  to  first  order  in  S/Rm, 


MRiClKlR^m  -  «,)  +  b/R^]  ~<  7*  >  BlI^„{SIR„)  +  r^a{dM fdr)dr  (5) 

Z  ''  Rm 


where  <  7,^  >=<  Bz  >  is  the  average  magnetic  pitch  angle  in  the  boundary  layer 
transition  zone  and  w,  =  Vlzj^K{Rm)-  GL  show  that  8  is  approximately  equal  to  the 
electromagnetic  screening  length.  We  will  not  make  this  assumption  unless  explicitly  noted. 
Except  for  the  second  term  on  the  left  hand  side  and  the  mass  loss  term  on  the  right  hand 
side,  eq.(5)  is  identical  to  eq.(43)  of  GL. 


Now  consider  the  mass  loss  term.  GL  treat  mass  flow  out  of  the  disk,  for  example 
by  vertical  mass  flow,  in  a  very  approximate  manner  by  introducing  a  “gate  function” 
g{r)  such  that  dM /dr  =  4:irpCzg{r),  with  g{r)  =  0  at  and  outside  the  outer  edge  of  the 
boundary  layer  and  g{r)  =  1  well  inside  the  boundary  layer,  c,  and  p  are  the  sound 
speed  and  density,  respectively.  For  simplicity  we  will  account  for  vertical,  mass  flow 
out  of  the  disk  by  introducing  a  parameter  such  that  the  radial  mass  flow  rate 
Mr(r)  =  (1  —  /0m)[r  —  Rm\MIS  -f-  in  the  transition  boundary  layer.  In  this  case  M  is 
the  radial  mass  flow  rate  in  the  disk  well  outside  the  transition  layer.  For  /5m  =  0,  there 
is  no  radial  mass  flow  through  r  =  Rm,  while  for  {3^  =  I  there  is  no  vertical  mass  loss. 
Evaluating  the  mass- loss  term  in  eq.(5)  in  this  manner  results  in 

MRiaK(R,.)[Hl  +  \ilRA  ~<  1*  >  (6) 

where  0  <  0m  <  ^  and  0m  =  {20m  +  l)/3.  Scaling  to  numerical  values  appropriate  to  QPO 
sources,  and  assuming  Bz  =  p/r^,  eq.(6)  gives 
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If  y5r„(l  —  u;,)  <<  \8I Rm  then  Rm  does  not  explicitly  depend  on  the  form  of  6jRm.  In 
this  limit  the  criterion  for  determining  Rm  is  the  same  as  that  given  by  Wang  (1987).  GL 
did  not  include  the  second  term  on  the  left  hand  side  of  eq.  (6).  For  values  of  w,  near  the 
critical  value  in  the  GL  theory  this  term  is  negligible  anyway.  However,  using  the  magnetic 
pitch  expression  of  Wang  to  evaluate  the  accretion  torque  (see  Appendix)  gives  critical 
values  of  ujg  much  closer  to  unity.  In  this  case  the  two  terms  on  the  left  of  eq.  (6)  can  be  of 
the  same  magnitude. 


3.2.  The  Magnetic  Pitch 

To  proceed  further  we  must  specify  <  7^  >  and  Rm.  First  consider  the  magnetic 
pitch.  Following  Wang  (1987),  we  assume  that  the  azimuthal  component  of  magnetic  field 
is  generated  as  a  result  of  the  relative  rotation  between  the  magnetosphere  and  the  disk 
and  is  limited  by  either  annihilation  of  toroidal  flux  near  the  midplane  of  the  disk  or  by 
expulsion  of  flux  through  the  effects  of  magnetic  buoyancy  (Parker  1955).  Setting  the  loss 
rate  equal  to  the  rate  of  generation  gives 

7^  =  B^IB,  ~  (A/C)'/2[(nK(r)  -  ^s)hlvA^\^'^  (8) 

where  vaz  =  Bzl{^'Kpyl'^  is  the  associated  Alfven  velocity.  (A/Q  is  a  dimensionless 
parameter  expected  to  be  greater  than  or  the  order  of  unity,  h  is  the  half-height  of  the  disk. 
Assuming  M  =  AitprVrh  in  the  disk  and  averaging  over  the  transition  zone,  eq.(8)  gives 

<H>=  =  7,  (9) 

where  u,  is  the  average  radial  flow  velocity  in  the  transition  zone.  In  the  original  treatment 
of  GL,  <  70  >  was  taken  to  be  a  constant  of  order  unity  in  the  transition  zone.  Also, 
to  calculate  the  aximuthal  magnetic  field  pitch  in  the  outer  disk  caused  by  interaction 
of  the  magnetosphere  with  the  Keplerian  flow  of  the  disk,  GL  estimated  the  size  of  B^ 
by  estimating  its  rate  of  amplification  by  differential  rotation  as  X\Uk  —  ^s\B^,  where 
A  is  a  numerical  factor  of  order  unity.  As  pointed  out  by  Wang  (1987)  and  Ghosh  and 
Lamb  (1991),  the  induction  equation  actually  implies  that  B^  is  generated  at  the  rate 
X\ilK  —  ^s\Bz.  This  expression  leads  to  eq.  (8),  which  generally  gives  smaller  values  of  the 
magnetic  pitch  in  the  outer  transition  zone  as  compared  with  GL  and  thus  gives  smaller 
values  of  the  spin-down  contribution  to  the  torque  and  therefore  larger  values  of  the  critical 
fastness  parameter.  In  any  case  70  is  likely  limited  to  values  ~  1  —  3  by  the  expansion  of 
the  magnetosphere  that  occurs  when  the  magnetospheric  magnetic  field  is  strongly  sheared 
(Ghosh  &  Lamb  1991).  From  the  numerical  estimate  given  below  we  find  that  <  70  >  given 
by  eq.(9),  for  parameters  appropriate  to  LMXB  Z  sources,  is  probably  below  this  range  of 
limiting  values. 
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3.3.  Width  of  Transition  Layer 

In  the  GL  theory  6  is  equal  to  the  electromagnetic  screening  length  given  by 

(c^/47r) 


8  = 


(VrCfeff) 


(10) 


where  Vr  is  the  average  radial  plasma  velocity  in  the  transition  layer  and  CTeff  is  the  effective 
disk  conductivity  for  poloidal  currents; 


C2  1 


CTeff  — 


7^ 


Combining  eq.(lO)  and  (11)  gives 


47r  hr  (^(r)  fl,) 

1  {VK\{h' 


(11) 


<7^  > 

Substitution  of  this  expression  into  into  eq.(7)  gives 

I'^r'ji.hl  Rm') 


(12) 


2/7 


R-  =  68.8km 


+  1(1/  <  70  >){VK/Vr){h/Rm)]- 


(13) 


In  their  treatment  GL  set  <  7^  >«  1  and  approximate  Vr  by  the  free-fall  velocity 
through  the  transition  layer; 

Vr  ~  {S/rY^‘^VK{l  -  =  Vff  (14) 

In  this  limit,  eq.(12)  and  (14)  give  [see  also  White  and  Stella  (1988)  and  GL  (1979)] 

{8lr)aL  —  {hlrYl^{\  —  u;4)^^^(l  (15) 


While  eq.(14)  is  ein  appropriate  estimate  when  w,  <<  1,  as  w,  — ♦  1,  the  characteristic 
radial  velocity  in  the  transition  layer  may  be  better  approximated  by  the  radial  drift  velocity 
in  the  disk  Vd  just  outside  the  transition  layer,  vj,  is  determined  primarily  by  the  viscosity 
of  the  disk.  For  a  radiation  pressure  dominated  disk  (RPD),  in  the  a-disk  approximation, 
(Treves,  Maraschi  &:  Abramowicz  1988) 

{vKlvd){hlr)  =  (1.4/ac)(2C?M/c2)-3/2(GM)^/V  (16) 

For  a  RPD  disk,  the  condition  that  the  radial  velocity  exceed  the  free-fall  value 
is  (1  -  uie)  «  5^^^5/2(^^/22.5km)-5/2.  Since  we  expect 

(1  —  a?,)  «  0.04  —  0.06  (see  below),  the  appropriate  limiting  value  for  Vr  depends  on  the 
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viscosity  parameter  a.  VVe  will  consider  this  further  below.  Assuming  for  the  moment 
that  Vr  =  Vd  and  conditions  appropriate  to  LMXB  Z-sources  {Bz  <  10^G(i?a/r)^  and 
(1  —  a?,)  «  0.06),  we  find  <  7^  >~  0.1a~^/‘‘(Bj/10^G)“^/^(M/10'®gm  Thus  it  is 

likely  that  the  average  magnetic  pitch  in  the  transition  zone  does  not  exceed  unity. 

Our  use  of  equations  describing  a  RPD  disk  implicitly  carries  the  assumption  that 
instabilities  do  not  destroy  the  disk  altogether.  There  is  good  reason  to  think  that  this 
is  the  case  when  HBO  are  present,  because  the  constancy  of  the  QPO  and  its  gradual, 
repeatable  variation  with  changing  luminosity  suggest  a  degree  of  constancy  associated 
with  the  disk  structure.  It  is  also  possible  that  instabilities  associated  with  RPD  disks  may 
lead  to  clumps  in  the  accretion  flow  which  are  a  necessary  ingredient  of  the  beat  frequency 
model,  since  if  there  is  no  azimuthal  variation  of  density  in  the  transition  zone  then  the 
QPO  will  be  absent. 

In  the  limit  where  the  transition  zone  velocity  is  the  free-fall  value,  eq.(13)  gives  the 
result  that  would  have  been  obtained  by  combining  the  GL  prescription  for  Rm  and  S/Rm 
with  the  RPD  disk  prescription  for  h/ Rm.  However,  eqn.(13)  is  more  general  than  this. 

In  order  to  explore  a  range  of  viable  magnetospheric  models  and  guided  by  the  above 
considerations  about  Vr  and  <  7,^  >,  we  adopt  the  following  parameterizations; 

<  >=  (Tw)*’  ,  V.  =  (17) 


where  7tu  is  given  by  eq.(9). 

The  plausible  range  for  the  exponents  and  is  0  </?-,<  1  and  0  <  <  1.  The 

GL  values  are  obtained  for  =  Q  and  =  0,  while  the  Wang  values  are  obtained  for 
=  1  and  =  1.  Eqs.  (9),  (12),  (14),  (15),  (16)  and  (17)  provide  a  consistent  set  of 
equations  for  determining  Rm  in  the  limit  of  a  RPD  disk  once  Bi(r)  is  specified.  We  assume 
that  Bz{r)  =  where  /ie//(r)  is  an  effective  magnetic  dipole  moment  that  includes 

corrections  due  to  screening  currents  in  the  disk. 

At  this  point  we  can  derive  a  relationship  between  the  spin  frequency  and  the  effective 
magnetic  dipole  moment,  utilizing  the  observational  constraint  of  the  beat  frequency  model 
implied  by  Eq.  (1). 


= « X 104 

(i/,  +  uqyi^  {2T:h[vz  +  vq)  27; 


Pv  1-/Jtp 

'mVd  ^ff 


(18) 


where  Mn  is  the  radial  mass  flow  outside  the  transition  region  in  units  of  lO^^gm  s 


-1 
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3.4.  Effective  Magnetic  Moment 

The  final  major  ingredient  in  our  model  is  a  treatment  of  the  modification  in  magnetic 
field  structure  caused  by  the  (partial)  exclusion  of  magnetic  flux  from  the  disk  and  the 
diamagnetic  behavior  of  the  neutron  star.  The  main  effect  is  an  enhancement  of  the 
magnetic  field  in  the  magnetosphere.  This  is  represented  by  an  effective  magnetic  moment: 

=  tx[Rml{Rm  -  Rs)]F  (19) 

where  Rg  is  the  equatorial  radius  of  the  neutron  star,  fj,  is  the  magnetic  moment  of  the  star 
in  the  absence  of  the  disk.  F  is  a  correction  factor  of  order  unity  that  characterizes  what 
fraction  of  the  field  lines  between  Rm  and  the  light  cylinder  are  excluded  from  the  disk. 

It  will  be  taken  to  be  unity  here.  Equation  (19)  is  the  simplest  functional  form  for  fi^ff 
that  has  the  correct  behavior  for  large  Rm  (i-e.  fXeff  and  also  conserves  magnetic  flux 

(when  F  =  1).  Equation  (19)  as  given  above  does  not  include  allowance  for  any  quadrupole 
terms  in  the  magnetic  field  of  the  bare  star  that  may  become  appreciable  within  a  few  tens 
of  kilometers  of  the  surface.  Such  a  quadrupole  term  can  increase  the  ratio  Hejjl to  a 
larger  value  than  that  given  by  (19). 


3.5.  Equilibrium  Spin  Period  Constraint 

The  prescription  for  the  magnetic  pitch  angle  given  by  eq.(8)  has  important 
consequences  for  the  accretion  torque  and  the  equilibrium  stellar  spin  frequency.  Wang 
(1987)  showed  that  for  a  gas-pressure  dominated  disk  the  accretion  torque  vanished  for  a 
critical  value  of  the  fastness  parameter  =  0.94  (compared  with  0.35  in  the  original  GL 
theory).  More  recently,  Wang  (1995)  found  for  various  physically  plausible  estimates  of  the 
dissipation  timescale  r,  that  a;,c  =  0.875  —  0.950.  A  similar  calculation  (See  Appendix)  can 
be  done  for  a  radiation  pressure  dominated  disk,  with  the  result  that  u;,c  Ci  0.918.  This 
result  should  be  regarded  as  a  lower  limit  since  in  a  system  containing  a  rapidly-spinning 
neutron  star  the  spin-down  contribution  to  the  torque  is  likely  cutoff  near  the  light- cylinder 
radius  in  the  disk,  Ric-  In  a  typical  LMXB  Z  source,  Ric  =  c/fij  may  only  be  a  few  times 
larger  than  the  magnetosphere  radius.  As  a  result  Uac  could  be  as  large  as  0.978. 

Within  the  context  of  the  BFMA  model,  the  stellar  spin  frequency  can  be  expressed  as 

Va  =  VQbJalil  —  UJs)  (20) 

For  the  LMXB  Z  sources,  such  as  GX5-1,  we  do  not  know  the  location  along  the  Z  diagram 
where  the  accretion  torque  is  zero.  If  we  assume  that  the  torque  is  zero  somewhere  along  the 
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normal  branch,  then  the  highest  observed  HB  QPO  frequency,  along  with  the  theoretical 
lower  limit  on  a;,c  discussed  above  (ie.,  0.918),  places  a  lower  limit  on  Ug.  For  GX5-1,  uq  as 
high  as  36  Hz  has  been  observed  on  the  HB.  This  implies  i/,  >  400//’^. 


4.  Comparison  With  QPO  Observations 


To  use  the  above  magnetospheric  model  to  fit  HBO  oscillations,  it  is  convenient  to 
assume  the  form 

i/Q  =  a{Lx)^  -  Ug  (21) 

with 

Lx  =  ^gravMacc  (22) 

^grav  is  the  gravitational  potential  at  the  surface  of  the  neutron  star  and  the  exponent  ^  is 
related  to  the  logarithmic  derivative  of  Rm  with  respect  to  M  by 


3  dlogRm 
2  dlogM 


(23) 


Mace  is  the  amount  of  mass  flow  which  ultimately  gets  accreted  onto  the  neutron  star.  In 
our  model,  this  is  related  to  M  and  ^rn-  If  =  1  then  there  is  no  vertical  mass  loss  and 
Macc  =  M.  If  the  corotation  radius  is  outside  the  transition  region,  i.e.  Rc  >  Rm  +  ^  then 
the  mass  flowing  vertically  out  of  the  disk  is  trapped  by  the  magnetic  field  and  is  eventually 
channeled  onto  the  neutron  star,  and  Macc  =  M  still  holds.  If  Rm  <  Rc  <  Rm  +  ^5  then  the 
vertical  mass  flow  outside  the  corotation  radius  is  ejected  from  the  system  by  the  propeller 
mechanism  ().  In  this  case  Macc  =  [(1  —  ~  l)i?m/<5  +  0m]M.  In  many  cases, 

simply  letting  Macc  =  M  is  sufficient.  However,  as  u;,  — »  1,  letting  Macc  =  ^mM  may  be  a 
better  approximation. 

Some  authors  at  this  point  introduce  an  additional  parameter  such  that  Lx  varies  as 
M^^.  For  simplicity,  we  refrain  from  introducing  such  a  parameter.  It  proves  unnecessary 
for  obtaining  satisfactory  agreement  with  observation. 

In  applying  this  form  to  HBO  data  it  is  worthwhile  to  recognize  that  there  is  an 
essentially  inverse  relationship  between  Ug  and  This  can  be  seen  in  a  simple  way  if  we 
assume  that  the  fit  is  constrained  to  pass  through  a  pair  of  points  on  the  curve  t/q  =  vq{Lx). 
Call  these  two  points  {vq,ii  Lx,i)i  i  =  1,2.  These  could  be  taken  to  be  representative  points 
near  ends  of  the  curve.  Then 


fiXu  -  1)  _ 


(24) 
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where  xl  =  Lx,2l Lx,\  and  Xi/  =  is  completely  determined  by  two  data  points 

and  if.  Given  f ,  the  allowed  range  of  for  a  source  is  determined  by  the  error  propagation 
in  the  above  equation.  Note  that  the  calibration  of  the  luminosity  does  not  enter  into 
determination  of  v^',  only  the  ratio  xl  affects  the  result.  Since  relative  uncertainties  in  the 
numerical  values  of  the  observational  quantities  in  the  equation  are  small  in  practice  for 
a  source  such  as  GX  5-1,  the  estimate  of  the  spin  frequency  in  such  a  source  is  largely 
dictated  by  the  estimate  of  f .  This  is  illustrated  in  Figure  2,  which  shows  allowable  regions 
in  I/,,  f  space  for  GX  5-1.  For  values  of  i/,  and  f  within  the  shaded  region,  an  acceptable  fit 
(which  we  define  as  z/q,  determined  from  eq.(21),  being  within  1  Hz  of  the  observed  vq)  of 
eq.(21)  to  the  data  for  GX5-1  can  be  obtained. 

f  can  be  derived  from  the  magnetospheric  model  presented  above.  The  general  result 
is  algebraically  complicated  and  not  directly  very  enlightening.  In  order  to  illustrate  the 
dependence  of  f  on  parcimeters  such  as  and  MRmi  in  Table  (2) 

we  present  analytic  expressions  for  f  for  two  sets  of  simple  limiting  cases. 

Since  we  expect  to  be  positive  [e.g.,  from  eq.  (19),  A^  =  (^e////z  —  1)  « 

(1  —  R,/Rm)~^  —  1],  the  enhancement  of  the  magnetic  field  within  the  magnetosphere 
by  screening  currents  flowing  within  the  disk  tends  to  “stiffen”  the  magnetosphere  (i.e. 
lower  the  value  of  f  for  a  wide  range  of  parameters.  For  low  Rm  responds  with  relative 
insensitivity  to  variations  in  mass  accretion  rate.  )  If  neutron  stars  in  LMXB  Z  sources  are 
rotating  as  fast  as  the  fastest  observed  radio  pulsars  (641  Hz),  then  A^  could  be  as  large  as 
2.4.  This  would  require  a  “stiff”  neutron  star  equation-of-state.  For  a  relatively  “soft”  eos, 
A^  ~  1.4. 

In  all  cases  except  (d)  in  Table  (2),  ^  either  does  not  depend  on  uz,  or  is  lowered  for 
increasing  values  of  Except  for  (c)  and  (e),  all  cases  allow  low  values  of  A  stiff 
magnetosphere  gives  good  fits  to  data  from  GX5-1,  for  comparatively  high  stellar  spin 
frequencies.  For  example,  with  a  stellar  spin  period  near  that  of  the  fastest  radio  millisecond 
pulsars  (1.5  ms),  the  GX5-1  HBO  requires  ^  ~  0.075. 

Given  a  value  of  ^  (i.e.,  an  assumed  spin  frequency),  the  QPO  observations  and  the 
magnetospheric  model  constrain  the  neutron  star  magnetic  moment  fi.  Table  (3)  shows  the 
values  of  /x  derived  for  the  limiting  cases  shown  in  Table  (2),  assuming  i/,  =  641  Hz  and 
(  =  0.075. 

Since  we  expect  (A/C)  >  1,  a  <  1,  and  A^  >  1,  we  conclude  that  the  model  yields 
magnetic  field  strengths  in  a  range  consistent  with  those  inferred  from  observed  spin-down 
rates  in  the  fastest  millisecond  radio  pulsars  (/X27  ~  0.1  ~  0-b)  (Camilo,  Thorsett  &  Kulkarni 
1994). 
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In  Figure  3  we  show  a  fit  (solid  line)  of  eq.  (21),  with  i/,  =  641  Hz  and  w  0.075,  to 
the  HB  QPO  of  GX5-1.  Also  shown  (dashed  line)  is  a  fit  of  eq.  (1)  to  the  data,  directly 
using  eq.(13)  to  determine  Rm  rather  than  resorting  to  the  power  law  (^  =  constant) 
approximation  used  above.  We  have  assumed  =  1,  /?„  =  1,  and  =  1.  Acceptable  fits 
can  be  obtained  for  a  relatively  broad  range  of  parameters  representing  cases  intermediate 
between  (a)  and  (e)  in  Table  (2). 

We  can  also  use  eqs.  (18)  and  (22)  to  directly  calculate  the  effective  magnetic  moment. 
For  the  case  where  jSm  =  l,/3y  =  1,  =  1,  we  find 


(25) 

where  Lx/e  =  Lxj Ledd  is  the  luminosity  in  units  of  the  eddington  luminosity.  If  we  let 
=  ^  while  keeping  /?„,  =  l,y9t,  =  1  then  we  obtain 


^^l}i27{Rm)  =  6.8  X  10^ 


x/e 


2.Qarn?^^  -f- 


500 


(26) 


Since  these  two  equations  are  for  the  effective  magnetic  dipole  moment,  we  expect 
the  intrinsic  dipole  moment  to  be  somewhat  less  than  this.  If  we  adopt  the  form  given 
by  eq.  (19),  we  expect  fi  ss  A^e//(1  —  6.7  x  10~^(i'j  +  (assuming  R,  =  10®cm).  In 

Figure  4  we  plot  eqs.  (25)  and  (26)  assuming  a  =  0.1,  (A/C)  =  1,  m  =  1.4,  i/q  =  36Hz  and 
Lx/e  =  1?  with  and  without  the  correction  for  magnetic  screening.  We  have  also  plotted 
several  known  magnetic  dipole  moments  of  millisecond  pulsars  versus  their  spin  frequency 
(Camilo,  Thorsett  &  Kulkarni  1994). 


5.  Discussion 

The  model  presented  here  represents  an  alternative  picture  of  the  HB  within  a  broader 
synthesis  of  Z  source  behavior  such  as  that  of  Lamb  (1989).  We  have  not  attempted  to 
treat  the  other  branches  of  the  Z  diagram  because  it  is  clear  that  a  different  mechanism  is 
responsible  for  QPO  oscillations  on  these  branches.  The  central  issue  under  investigation 
has  been  the  relationship  between  i/,  and  C,  represented  in  Figure  2,  and  the  constraints 
that  HB  QPO  place  on  these  parameters.  We  note  the  following  considerations  favoring 
our  modified  beat  frequency  model:  (i)  incorporation  of  relevant  physical  effects  such  as 
radiation  pressure;  (ii)  consistent  treatment  of  the  effects  of  stellar  spin  in  determining  the 
accretion  torque  and  the  magnetosphere  radius;  (iii)  preservation  of  the  fit  of  the  model 
to  the  observed  luminosity  dependent  QPO  frequency;  and  (iv)  agreement  of  stellar  spin 
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periods  and  magnetic  moments  determined  by  fits  to  the  model  with  those  observed  for  the 
fastest  radio  pulsars.  The  central  conclusion  we  reach  is  that  U3  in  Z  sources  is  at  least 
~  400  Hz  and  that  values  as  high  as  observed  in  the  fastest  millisecond  radio  pulsars  (641 
Hz)  are  consistent  with  observation  and  magnetosphere  theory. 

Recent  observations  of  Sco  X-1  and  (van  der  Klis  et  al.  1996)  4U1728-34  (Strohmayer, 
Zhang,  and  Swank  1996)  with  XTE  have  revealed  ~  1000  Hz  and  ~  800  Hz  QPO 
respectively,  which  may  represent  the  Keplerian  frequency  of  the  inner  edge  of  the  disk. 
Although  Sco  X-1  has  not  been  observed  on  the  Horizontal  Branch,  and  therefore  has  no 
measured  HB  QPO,  if  we  assume  that  u,  ~  a;sc  ~  1  then  the  spin  frequency  would  also  be 
~  1000  Hz.  Noting  that  -j-  vq  in  our  model,  we  can  estimate  the  surface  magnetic 

field  H  ~  1  —  3  X  10®  G,  consistent  with  the  millisecond  pulsars.  The  question  which  still 
has  to  be  answered  is  why  the  coherent  pulse  period  has  not  yet  been  seen,  while  the  inner 
edge  of  the  disk  has  been  seen.  One  possible  explanation  is  that  a  weak  pulse  from  the 
neutron  star  surface  would  be  strongly  attenuated  by  gravitational  lensing  as  discussed 
in  section  2,  while  the  QPO  from  the  disk  at  a  distance  of  several  stellar  radii  would  be 
essentially  unaffected.  This  weak  pulsar  signal  could  then  be  obscured  by  the  broad  (~  50 
Hz  FWHM)  disk  QPO.  Further  analysis  of  the  data  may  yet  reveal  an  underlying  coherent 
signal. 

In  a  subsequent  paper,  we  will  consider  further  the  robustness  of  the  foregoing 
model.  For  the  present,  we  will  simply  mention  some  of  the  observational  and  theoretical 
considerations  that  support  the  revised  model  which  will  be  discussed  in  detail  later. 
First,  the  width  of  the  QPO  features  and  their  variation  of  intensity  is  a  natural  result  of 
the  model.  Observations  of  HB  QPO  from  GX  17-f-2  have  suggested  that  the  dominant 
mechanism  contributing  to  the  width  of  the  QPO  is  smearing  in  uq  due  to  repeated 
interaction  of  accreting  clumps  with  stellar  magnetic  field  as  the  clumps  drift  inward  across 
a  transition  zone  of  finite  radial  extent  (Penninx  et  al.  1990).  The  above  expression  for  Rm 
implies  a  spectral  width  more  consistent  with  observations  of  GX5-1  (van  der  Klis  et  al. 
1985)  than  that  derived  from  the  conventional  BFMA  model.  Secondly,  the  relationship 
between  QPO  frequency  and  the  evolution  of  the  neutron  star’s  magnetic  field  derived  from 
the  model  is  consistent  with  the  observed  dispersion  in  the  magnetic  fields  inferred  in  the 
fastest  millisecond  radio  pulsars  and  the  6  <  uq  >(<  uq  >  observed  in  Z-sources. 

A  possible  argument  against  the  model  is  the  high  shot  rate  required  for  the  implicit 
correlation  between  the  QPO  and  the  low  frequency  noise  component  of  the  Fourier  power 
spectrum.  Lack  of  detection  of  such  a  correlation  in  GX5-1  requires  a  shot  rate  greater  than 
400Hz,  which  has  been  claimed  to  be  inconsistent  with  the  BFMA  model  (Norris  et  al.  1990). 
However,  that  calculation,  which  implied  a  maximium  shot  rate  of  300  Hz,  assumed  no 
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variance  in  the  shot  amplitude.  However,  if  there  is  a  distribution  of  shot  amplitudes,  such 
as  an  exponential  distribution,  the  limit  of  the  shot  rate  becomes  600  Hz,  consistent  with 
observation.  Another  possible  problem  is  the  onset  of  the  Chandrasekhar- Friedman- Schutz 
(CFS)  gravitational  radiation  reaction  instabilility  when  the  neutron  star  is  spun  up  to  the 
frequencies  (i'j  >  640Hz)  implied  by  the  model.  However,  a  calculation  by  Meidell  and 
Lindblom  (1991)  suggests  that  if  the  neutron  star  is  below  10®  K  and  in  a  superconducting 
state,  then  the  viscosity  due  to  electron  scattering  from  superfluid  vortices  will  suppress  the 
CFS  instability. 


6.  SUMMARY 

We  have  presented  a  new  model  of  horizontal  branch  QPOs  that  is  based  on  the 
beat  frequency  modulated  accretion  model.  Its  departures  from  the  original  BFMA  model 
correct  recognized  deficiencies.  The  model  predicts  high  stellar  spin  frequencies  (  >  400  Hz) 
that  fall  in  a  range  where  pulsation  upper  limits  pose  no  problem.  The  magnetic  field  and 
spin  frequency  range  are  also  compatible  with  those  of  the  fastest  radio  pulsars.  However, 
detection  of  the  stellar  spin  period  is  the  only  hope  of  rigorous  confirmation.  Short  of  that, 
the  model  presented  comes  close  to  incorporating  all  available  information. 
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A.  Accretion  Torque  and  the  Critical  Fastness  Parameter 

Assuming  the  magnetic  field  pitch  distribution  7^  given  by  eq.(8)  and  a  thin, 
gas-  pressure  dominated  disk,  Wang  (1987)  calculated  the  accretion  torque  and  found 
Uac  =  0.941,  where  iOsc  is  the  critical  fastness  parameter  for  which  the  accretion  torque 
vanishes  (Ghosh  &  Lamb  1979).  In  the  GL  model,  u;,c  =  0.35. 

Observations  of  QPO  in  the  x-ray  pulsars  Cen  X-3  (Tennant  1988,  Takeshima  et  al. 
1991),  EXO  20304-375  (Angelini,  Stella  and  Parmar  1989)  and  X1627-673  (Shinoda  et  al. 
1990)  provide  important  observational  tests  of  both  the  beat  frequency  model  of  QPOs 
and  of  accretion  torque  models  (Shibazaki  1989).  We  note  again  that  observations  in  these 
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systems  are  not  necessarily  directly  relevant  to  the  QPO  observed  in  LMXBs,  but,  at  least 
provide  supporting  evidence  for  some  of  the  key  assumptions  made  in  modelling  the  HB 
QPOs  in  the  LMXB  Z  sources. 

Table  (Al)  lists  the  observed  spin  frequencies,  QPO  frequencies  and  spin  period 
derivatives  for  each  of  the  above  sources.  Also  listed  are  the  values  of  the  fastness  parameter 
derived  by  assuming  the  beat  frequency  relation  of  eq.(l)  is  valid.  For  Cen  X-3  and 
X1627-673  the  fastness  parameter  is  much  larger  than  the  critical  fastness  of  GL  but  below 
that  of  Wang.  Both  of  these  sources  are  observed  to  be  spinning-up,  consistent  with  the 
Wang  model  but  not  with  the  GL  model.  For  EXO  2030+375,  the  fastness  parameter  is 
below  the  critical  value  for  both  the  GL  model  and  the  Wang  model. 


In  the  Wang  theory,  the  period  derivative  is  given  by 

w.  \ /1.4M0\473/io6cm\  , 

-p/p  =  4.28  X  io-VK)y  (i035i;^)(-^)  (-«-)>'''  (Al) 


/(a,.)  =  u.'/^  + 


I  —LO.— 


U> 


,3/2 


(A2) 


For  parameters  derived  from  observation  and  the  beat  frequency  relation,  eq(Al)  can 
be  used  to  calculate  (P/P),  which  can  then  be  compared  with  the  observed  value.  The 
theoretical  values  are  shown  in  the  last  column  of  Table  (Al).  Given  the  uncertainties 
in  determining  the  luminosity  (distance)  and  uncertainties  in  the  masses  and  radii  of  the 
neutron  stars  in  these  systems,  the  agreement  of  observation  with  theory  is  satisfactory.  In 
the  case  of  Cen  X-3  the  theoretical  value  of  P  differs  from  observation  by  a  substantial 
amount.  However  in  this  source  P  shows  considerable  fluctuations  about  the  mean  value 
shown  in  the  Table  (Al).  We  also  note  that  Parmar  et  al.  (1989)  analyzed  an  EXOSAT 
observation  of  EXO  2030+375  during  an  x-ray  outburst  in  which  the  luminosity  changed 
by  more  than  10^.  They  found  that  -P  oc  xhis  is  also  in  reasonable  agreement 

with  the  luminosity  dependence  expected  from  eq.(Al),  namely  —P  oc  The  GL  theory 
predicts  -P  oc  Similarly,  Wilson  et  al.  (1993)  have  found  that  for  GS0834-430 

—P  oc 

In  the  remainder  of  this  appendix  we  evaluate  the  accretion  torque  and  o/jc  in  the  limit 
of  a  radiation  pressure  dominated  inner  disk.  Following  GL  and  Wang,  the  accretion  torque 
is  given  by 

N  =  \B^\B,r^dr  -  \B4,\B,r^dr  (A3) 

jRm  JHc 

where  Rc  =  (GAf/f)^)'/^  is  the  corotation  radius  and  Ro  is  either  the  outer  radius  of  the 
disk  or  ~  the  light- cylinder  radius.  Magnetic  field  penetrating  the  disk  for  r  >  Rc  provides 
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a  spin-down  torque,  while  field  between  Rm  and  Rc  enhances  the  spin-up  torque.  As  first 
shown  by  GL,  the  spin-down  contribution  to  the  torque  can  equal  or  exceed  the  spin-up 
contributions. 

The  spin-down  contribution  to  N  can  be  written  as 

N-  =  —  [  'y^{r)Blr^dr  (A4) 

JRc 

or,  using  eq.(8), 

iV_  =  [{r/Rcff'^ -l]^^^[{M/r){h/r){vK/vr)^lKY^‘*Bl^'^r^dr  (A5) 

For  a  RPD  disk  we  assume 

{h/r){vK/vr)  =  {\Alac)Rfl'^{GMfl'^r  (A6) 

Substituting  eq.(A6)  into  (A5)  and  assuming  Bz  =  /z/r^, 

N.  =  {y  -  (A7) 

where 

Nc  =  (A8) 

The  integration  variable  y  =  {rj Rc)^!^  and  yo  =  {Roj Rc)^^^ .  The  spin-up  torque  can  be 
written  in  a  similar  form; 

iV+  =  M{GMRmY^^  +  \Nz  Ai  -  yY'h'^^'^'^dy  (A9) 

O  Ju)s 

noting  that  w,  =  {Rm! Ref  1"^. 

Now  define  =  M^GMRm-  It  can  be  shown  that 

=  47V2(1  -  u;z)u:f/\6/R„,f[8/Rm  +  2fm{l  -  u;,)]'"  (AlO) 

Following  GL,  it  is  convenient  to  introduce  the  dimensionless  accretion  torque  defined 

by 

N  =  Neviu^f  (All) 

If  we  approximate  as  y~^,  then  the  integrals  in  eq.(A7)  and  eq.(A9)  can  be  done 

analytically  with  the  result  that 

(taa-‘(!>o  -  1)'^^  -  ‘^)) 

(A12) 


,(u..)  =  l+i(l  +  r(l-c.<,)''=|[u.y-+ 


a;®/‘*(tanh  ^(1— 


(1 
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where  T  =  2$rn.{iK)''^^[{hl Rm){‘>JKlvr)]~^^\{M ! Note  that  T  does  aot 
explicitly  depend  on  the  stellar  spin  frequency.  It  does  depend  on  Rm- 

The  Wang  magnetosphere  radius  criteria  corresponds  to  the  limit  F  =  0.  To  determine 
the  critical  fastness  parameter  oJac  for  which  ?/  =  0,  Wang  also  assumed  2/0  =  00  and  a  gas 
pressure  dominated  disk.  He  found  u^c  =  0.94.  Evaluating  eq.(Ar2)  for  F  =  0  and  yo  =  00 
gives 

riM  =  1  -  (t/6)(1  -u,,)-'/^tanh-'((l  -  u-.)'/")  (A13) 

This  function  is  shown  as  the  F  =  0  curve  in  Figure  (Al).  For  w,  =  0.918,  the  accretion 
torque  vanishes.  This  value  is  slightly  less  than  the  critical  fastness  for  a  gas  pressure 
dominated  disk.  The  other  curves  in  Fig.(Al)  show  the  behavior  of  t){u!s)  for  fixed  values 
of  F.  For  parameters  appropriate  to  RPD  disks  and  typical  LMXB  Z  sources,  we  expect 
F  <  10.  As  F  increases  from  zero  the  critical  fastness  parameter  decreases  but  always 
satisfies  the  condition  u;,c  >  0.65. 

For  a  very  rapidly  spinning  neutron  star,  the  light  cylinder  radius, 

Ri^  =  c/Cla  =  74.45km(640Hz/i/4),  may  only  be  a  few  times  the  corotation  radius 
Rc  since 

Ric/Rc  =  =  3.2(640Hz/i/,)2/3(i.5M©/M)‘/^  (A14) 

In  this  case  letting  yo  00  is  not  the  appropriate  limit  and,  in  any  case,  the  integrand  in 
eq.(A7)  becomes  a  poor  approximation  when  r  approaches  or  goes  beyond  Ric-  Lacking 
a  theory  valid  in  this  regime,  we  resort  to  approximating  the  effects  of  rapid  rotation  by 
taking  the  upper  limit  of  integration  in  eq.(A12)  as 

Vo  =  {Rid  Ref =  ud'^GMIR^dy^l^  (A15) 

The  curves  in  Figure  (A2)  show  r]{ijjs)  for  various  values  of  F  assuming  relation  (A14). 
The  main  effect  of  introducing  this  cutoff  is  to  increase  uJsc-  For  F  =  0,  uJsc  =  0.978.  As  F 
increases,  decreases  but  satisfys  the  inequality  >  0.88. 

Figure  (A3)  is  a  contour  plot  illustrating  the  dependence  of  rj  on  a;,  and  F  for  two 
limiting  cases;  (a)  =  oo  and  (b)  yo  =  {Rid In  both  cases  ?/  =  1  for  a;,  =  0.  ^ 


^Recall  that  we  approximated  by  y^.  This  approximation  breaks  down  for  w,  =  0.  If  this 

approximation  is  not  made,  the  value  of  y  is  slightly  greater  than  1  for  Us  =  0.  In  evaluating  the  torque  in 
the  gas  pressure  dominated  case,  Wang  (1987)  makes  a  similar  approximation;  y29i/i20  ^  ^^5/2 
the  approximation  leads  to  a  divergence  of  rj  for  =  0.  If  the  approximation  is  not  made  then  t)  is  well 
behaved  for  =  0.  Ghosh  and  Lamb  (1991)  have  criticized  Wang^s  torque  theory  because  of  the  presence 
of  this  divergence.  This  criticism  is  not  valid. 
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Table  1.  Limits  on  pulse  fractions  in  Z  sources  for  spin  frequencies  up  to  500  Hz.“ 


Source  pulse  fraction 


Sco  X-1 

<  0.57% 

GX  340+0 

<  1.8% 

GX  17+2 

<  1.6% 

Cyg  X-2 

<  1.4% 

GX  5-1 

<  0.99% 

GX  349+2 

<  1.3% 

^99%  confidence  Vaughan 
et  al.  1994. 


Table  2.  The  exponent  ^  for  several  limiting  cases  indicated  by  the  values  of  the 

parameters  in  the  first  column. 


parameter  values 


range  of  ^ 


(a)  =  1; /?v  =  1  3/7(1  + 4/7A^  +  2/7u;,/(l-w,))-i  0  <  ^  <  0.102 

(b)  /?^  =  1;  /?„  =  0  3/10(1  +  8/15A;.  +  3/15w,/(l  -  u>,)  +  l/20w,/(l  +  0  <  ^  <  0.092 

(c)  =  0; /?„  =  0, 1  3/7(1 +  4/7 A;,)- 1  0.429  >  ^  >  0.181 

0^{l-U,)»l/26/Rm 

(d)  =  1, 0;  /?„  =  0  3/25(1  +  U/25\^  -  6/25w,/(l  -  w,)  +  3/25w,/(l  + 

(e) ^^  =  l,0:/?„  =  1  3/5(1 +  4/5A^)-‘ 


^  >  -0.071 
0.205  <  ^  <  0.60 
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Table  3.  Neutron  star  magnetic  moment  for  various  cases  defined  in  Table  (1),  assuming 
'  that  u,  =  641  Hz  and  ^  =  0.075. 


Case 

(a) 

6.35(A/C)-^/^a^/®(A^*  + 1)”^ 

(b) 

0.744(A/C)-®/*(A^  + 1)-' 

(d) 

0.233;5V«(A^  +  l)-i 
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Table  Al.  QPO  properties  of  x-ray  pulsars.  Shown  in  the  table  are  the  observed  spin  and 
QPO  frequencies,  the  fastness  parameter  (assuming  the  beat  frequency  relation),  the 
observed  spin  period  derivative,  and  the  theoretical  spin  period  derivative  calculated  using 
eq.(Al).  Parameters  in  the  first  three  columns  are  summarized  in  Shibazaki  (1989). 


Pulsar 

i/s  (Hz) 

UQ  (Hz) 

^3 

(P/FU(yr-‘) 

{P/P)theoTy{yr  ^) 

SMC  X-1 

1.410 

0.010 

0.993 

-f2.2  -  +2.6  X  10-3 

4U0115-I-63 

0.277 

0.062 

0.817 

-4.3  X  lO-'* 

V0332-K53 

0.229 

0.051 

0.818 

-5.7  X  10-“ 

Cen  X-3^ 

0.207 

0.035 

0.855 

-3  X  10"* 

-1.9  X  10-3- 

-0.6  X  10-“ 

X1627-673'’ 

0.131 

0.04 

0.76 

-2.2  X  10-^ 

-4.8  X  10-“ 

EXO  2030-^375' 
A0535-b26 

0.0239 

0.0096 

0.21 

0.10 

-3  X  10-2 

-6  X  10-2 

*The  maximum  (minimum)  x-ray  luminosity  is  8  x  10^^(2.6  x  10^®)erg  s~*  for  an  assumed 
distance  of  10  kpc  (Bradt  &  McClintock  1983).  {PfP)obs  shows  considerable  fluctuations 
with  time  about  the  value  shown  (Nagase  1989).  The  range  shown  for  {PIP)theory 
corresponds  to  the  range  of  observed  luminosity. 

“^The  x-ray  luminosity  is  lO^’^erg  s~^  for  an  assumed  distance  of  8  kpc  (Shinoda  et  al. 
1990).  (P/P)o6»  is  relatively  steady  in  this  source  (Nagase  1989). 

*^(P/P)o6,  is  from  Parmar  et  al.  (1989)  for  a  luminosity  of  9  x  lO^’^erg  s~^,  assuming  a 
distance  of  5  kpc. 

In  all  cases  the  neutron  stax  mass  and  radius  are  assumed  to  be  IAMq  and  10®  cm, 
respectively. 
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Fig.  1. —  Comparison  of  pulsation  upper  limits  for  various  LMXB  sources  (solid  curves 
labeled  Sco  X-1  and  GX5-1;  other  LMXB  limits  mentioned  in  the  text  lie  between  these  two 
sources)  with  attenuation  expected  from  scattering  in  a  hot  corona  surrounding  the  source 
(dashed  curves).  The  coronal  attenuation  depends  on  Rc  and  r.  At  low  frequencies  the 
attenuation  (the  parameter  Rct)  is  constrained  by  the  observed  QPO  power  (cross-hatched 
region).  At  high  frequencies  the  attenuation  approaches  .  From  observed  temporal  lags 
and  spectral  fits  to  coronal  models,  the  Thompson  depth  is  likely  to  be  in  the  range  from  4 
to  7. 

Fig.  2. —  Allowable  regions  in  i/,,  ^  space  that  give  acceptable  fits  to  the  frequency  dependent 
QPO  observed  in  GX  5-1. 

Fig.  3. —  Comparison  of  a  fit  (solid  curve)  of  eq,(21),  with  i/,  =  641  Hz  and  (  =  0.075,  to 
the  HB  QPO  observed  in  GX  5-1.  The  dashed  curve  is  a  fit  of  eq.(l)  to  the  data,  directly 
using  the  magnetosphere  radius  relation,  eq.(16),  with  =  !,/?„  =  1,  and  A^^  =  1. 

Fig.  4. —  Magnetic  dipole  moment  as  a  function  of  spin  frequency  from  eqs.  (25)  (solid 
curve)  and  (26)  (dashed  curve)  assuming  a  QPO  frequency  of  36  Hz,  viscosity  parameter 
a  =  0.1,  neutron  star  mass  of  IAMq,  and  Eddington  luminosity.  If  we  correct  for  magnetic 
screening,  we  obtain  the  short  dashed  (25)  and  dotted  (26)  curves.  Also  plotted  are  the 
intrinsic  magnetic  moments  for  the  fastest  millisecond  pulsars  (dots)(Camilo,  Thorsett  & 
Kulkarni  1994). 
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Figure  (A3).  Contour  plots  of  diraensionless  toit)ue  ij,  illustrating  the  dependence  on  Fand  (0^. 
For  case  (a),  yo  =  «»,  while  for  case  (b),  y©  = 


